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Abstract 

' An integrable theory is developed for the perturbation equations engendered from 

, small disturbances of solutions. It includes various integrable properties of the per- 

turbation equations: hereditary recursion operators, master symmetries, linear repre- 
sentations (Lax and zero curvature representations) and Hamiltonian structures etc. 
and provides us a method to generate hereditary operators, Hamiltonian operators 
. and symplectic operators starting from the known ones. The resulting perturbation 

equations give rise to a sort of integrable coupling of soliton equations. Two examples 
Tj" , (MKdV hierarchy and KP equation) are carefully carried out. 

o 

1 Introduction 

Integrable nonlinear wave or evolution equations (for instance, KdV, NLS, SG and KP 
equations) are ideal mathematical models of real physical phenomena although they play 
an outstanding role in physical problems. Therefore for these equations we often need to 
take into account the effect of small perturbation so that their applicability may be ex- 
tended to higher order nonlinearity or larger amplitude waves. There are mainly two kinds 
of perturbation worthy studying for integrable equations. The one is the perturbation sit- 
uation of integrable equations themselves and the other one, the perturbation situation of 
solutions of the original unperturbed integrable equations. They all can provide approxi- 
mate solutions to real physical problems. 

In the context of soliton perturbation, what one considers usually is the first kind of 
perturbation situation. There have been quiet a few of powerful techniques for dealing with 
this kind of the perturbation (see Q || Q [J|| || and references therein). Various per- 
turbed cases of integrable equations have been considered, including the perturbed KdV 
and MKdV equations® H @ , the perturbed nonlinear Schrodinger equation® H , the per- 



turbed Burgers equation!!^! and the perturbed Benjamin-Ono equation'^:' etc. A detailed 
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survey for this kind of perturbation theory based upon the inverse scattering transforma- 
tion was provided by Kivshar and Malomed [O] . However there don't exist so many works 
devoted to the second kind of perturbation situation. Among them there are the following 
several works. Tamizhmani and Lakshmanan considered the complete integrability of the 
perturbed equations by small disturbance of solutions of KdV equation^ . We analyzed 
the similar perturbation situation of the whole KdV integrable hierarchy and pointed out 



a mistake in Ref. |T^] for the Hamiltonian structure E-i . Recently it has been observed 
by Kraenkel et al. that a multiple time scale expansion may relate the solutions of long 



surface water-waves and the Boussinesq equation to KdV integrable hierarchyl!3 lis . 

In this paper we would like to consider the second kind of perturbation. Mathemati- 
cally, this kind of perturbation yields interesting results. For example, we shall show that 
it preserves complete integrability. In other words, the equations generated by the per- 
turbation are still integrable and thus give rise to new examples of integrable equations. 
Moreover they are all special integrable coupling of the original integrable equations. 

We now introduce our notation and conception. Some notation comes from Refs. [jl7| , 



[ [Tql , [19]. Let M = M(u) be a suitable manifold possessing a manifold variable u (we 
write u as a column vector) and Mn = Mn{t)n) be another suitable manifold possessing 
a manifold variable fjiy = (t]q, ryf, • • • , T)n) t , N > 1, where rji, < i < N, are column 
vectors and T means the transpose of matrices. Assume that T(M),T(Mn) denote the 
tangent bundles on M and Mat, T*(M),T*(M^) denote the cotangent bundles on M 
and Mn, and C°°(M), C°°(Mn) denote the spaces of smooth functions on M and Mn, 
respectively. Throughout this paper we require that the column vector rji (0 < i < N) has 
the same dimension as the column vector u. Therefore we have r\i G R q (0 < i < N) when 
u G R q . Further let TJ(M) be the s-times co- and r-times contravariant tensor bundle and 
(TJ) U (M), the space of s-times co- and r-times contravariant tensors at u € M. We use 
X(u) (not X\ u ) to denote a tensor of J £ TJ(M) at u G M but sometimes we omit the 
point u for convenience while there is no confusion of the symbols. Note that linear maps 
$ : T(M) -> T(M), * : T*(M) -» T*(M), J : T*(M) -> T(M), G : T(M) T*(M) 
may be identified with the second degree tensor fields 2$ G T}(M), X* G T}(M), Tj G 



T 2 (M), T e G T 2 °(M) by the following relations &\ 

T$(u)(a(u),K(u)) =< a{u),<&{u)K{u) >, a G T*(M), K G T(M), 
Tq,(u)(a{u),K(u)) =< ^(u)a(u),K(u) >, a G T*(M), K G T(M), 
Tj(u)(a{u),p{u)) =< a{u),J(u)/3(u) >, a,/3 G T*(M), 
T e (u)(K(u),S(u)) =< e(u)K(u),S(u) >, K, S G T(M), 

where < •, • > denotes the duality between cotangent vectors and tangent vectors. 

A basic tool to handle various kinds of tensor fields is the conception of the Gateaux 
derivative. For a tensor fields X G TJ(M), its Gateaux derivative at the direction S G 



2 



T(M) is defined by 

X>(U)[S(U)]=?- X{U + £S{U)) 



(1.1) 

e=0 



de 

For four kinds of operators between the tangent bundle and the cotangent bundle, their 
Gateaux derivatives may be given similarly or by means of their tensor fields. The com- 
mutator of two vector fields K,S G T(M) and the adjoint map adx are defined by 

[K,S](u) = K'(u)[S(u)} - S'(u)[K(u)], ad K S = [K,S]. (1.2) 

The conjugate operator of an operator between the tangent bundle and the cotangent 
bundle is established in terms of the duality between cotangent vectors and tangent vectors. 
For example, we may calculate the conjugate operator <J>t : T*(M) — ► T*(M) of an operator 
$ : T(M) -> T(M) through 

< &(u)a(u),K(u) >=< a(u),<f>{u)K(u) >, a G T*(M), K G T(M). 

If an operator J : T*(M) -> T(M) (or : T(M) -> T*(M)) plus its conjugate operator 
equals to zero, then it is called skew-symmetric. 

Definition 1.1 For # G C°°(M) ; its variational derivative ^ G T*(M) is defined by 

< ^-(u),K(u) >=< 6 J^,K{u) >= H'{u)[K(u)}, K G T(M). 
ou ou 

If for 7 G T*{M) there exists H G C°°(M) so that 

H'{u)[K(u)] =< j(u),K(u) >, for all K G T(M), 
holds, then 7 G T*(M) is called a gradient field with a potential H G C°°(M). 

A cotangent vector field 7 G T*(M) is a gradient field iff 



(d 7 )(n)(K(n),S(n)) 

:=< i(u)[K(u)\, S(u) > - < Y(u)[S(u)],K(u) >= 0, K, S G T(M). 
If 7 G T*(M) is gradient, then its potential F G C°°(M) is given by 

H(u) = [ < j(\u),u > d\. 
Jo 



(1.3) 



Definition 1.2 ^4 linear operator $ : T(M) — > T(M) is called a recursion operator of 
Ut = K(u), K G T(M), if we have for all S G T(M) and u £ M 

*g>s<«) + s'M^SM - if («)[.(.)«(«)] + *(«)if'W[SM] - 0. (i.4) 
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Evidently a recursion operator : T(M) — > T(M) of it^ = if (u) maps symmetries into 
new symmetries of itt = -fT('u). 



Definition 1.3 A linear operator $ : T(M) — > T(M) is called a hereditary operator $1 
if the following equality holds 

<S>'(u)Mu)K(u)]S(u) - <f>(u)&(u)[K(u)}S(u) 

-<S>'(u)[$(u)S(u)]K(u) + $(u)&(u)[S(u)]K(u) = { ' ' 

for all vector fields K, S G T(M). 

When an evolution equation u t = K(u) possesses a time-independent hereditary recursion 
operator <I>, a hierarchy of vector fields <& n K, n > 0, are all symmetries and commute with 
each other. If the conjugate operator \P = $t of the hereditary operator : T(M) — ► T(M) 
maps a gradient field 7 G T*(M) into another gradient field , then \E' n 7, n > 0, are all 
gradient neldsOH. 

Definition 1.4 ^4 linear skew- symmetric operator J : T*(M) — > T(M) is called a Hamil- 
tonian operator if for all a,/3,^ G T*{M) we have 

< K(u), J'{u)[J(u)S(u)]T(u) > +cycle(K, S, T) = 0. (1.6) 

The corresponding Poisson bracket is defined by 

{H x ,H 2 }j{u) =< 6 -^l(u),J(u) 6 -^(u) >, H U H 2 G C°°(M). (1.7) 
ou ou 

A pair of operators J,M : T*(M) — > T{M) is called a Hamiltonian pair if J + cM is 
always Hamiltonian for any constant c. 



When J : T*(M) -> T(M) is Hamiltonian, we have^S 

jl-{H 1 ,H 2 }j=[J^,J^}, Hi, H 2 G C°°(M). 

Moreover if J, M : T*(M) — ► T(M) is a Hamiltonian pair and J is invertible, then <J> = 
Mr 1 : T(M) -> T(M) is hereditary EZ I . 

Definition 1.5 ^4 linear skew- symmetric operator : T{M) — > T*(M) is called a sym- 
plectic operator if for all K,S,T G T(M) we have 

< K(u),Q'(u)[S(u)]T(u) > +cycle(K,S,T) =0. (1.8) 

If : T(M) — » T*(M) is a symplectic operator, then its second degree tensor field Te G 
T2 (M) may be expressed as 

T e = d7 with < -y(u),K(u) >= f 1 < Q{Xu)Xu,K{u) > dX, K G T(M), 

Jo 

where dj is defined by (|L^). It is not difficult to prove that the inverse of symplectic 
operators are Hamiltonian if they exist and vice verse. 
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Definition 1.6 An evolution equation Ut = K(u), K S T(M), is called a Hamiltonian 
equation if there exists a function H G C°°(M) so that 

()H 

u t = K{u) = J{u) — {u). (1.9) 
du 

It is called a hi- Hamiltonian equation if there exist two functions H±,H2 £ C°°(M) and a 
Hamiltonian pair J,M :T* (M) — > T(M) so that 

u t = K(u) = J{u)°—±{u) = M°—^[u). (1.10) 
ou du 

There is another kind of Hamiltonian equations, which may be defined by symplectic 
operators. However, the above definition is more advantageous. For a bi-Hamiltonian 
equation above, there exist several beautiful characteristics in the aspects of algebra and 



geometryEl Jf** 

In this paper, we shall analyze the perturbation equations of the evolution equation 

u t = K(u), K G T(M) (1.11) 

under the perturbation series 

N 

un = Y. £ V VN = {Vo, Vi----, vlf, N > 1 (1.12) 

and their integrable properties. The paper is organized as follows. In Section 2 we pro- 
pose three useful theorems to generate Hereditary operators, Hamiltonian operators and 
symplectic operators from the perturbation series ( l-12j ) in a natural and explicit way. We 



show in Section 3 that the perturbation ( 1,12; ) preserves complete integrability, i.e. we want 



to show that the perturbation equations are still integrable equations provided that the 
equation under consideration is integrable. In Section 4, we apply the resulting theory to 
MKdV hierarchy and KP equation as illustrative examples. Section 5 contains some con- 
cluding remarks, where we give another perturbation series and compare its corresponding 



results with ones by (1.12) 



2 Hereditary, Hamiltonian or symplectic operators by per- 
turbation 

We make a perturbation series for any m > 

m 

where rji, < i < m, are all column vectors possessing the same dimension as u. We first 
analyze a few of properties of tensor fields and then establish three useful theorems to 
construct hereditary, Hamiltonian or symplectic operators in terms of a perturbation. 
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Lemma 2.1 We have for any X G TJ(M) 

d k X(ui) d k X(i 



—, i,j >k>0. 



de k de k 
Proof: Let j > i without loss of generality. Then we have 



(2.2) 



Ua = U, 



=i+l 



Further we get 



X(u J )=X(u i ) + o{e i ), 
from which the required equality follows. # 
Let us assume for any X G TJ(M) 

1 &X(uj) 



(*(«))% 



X« = =(X(u))«(r),), j > i > 0. (2.3) 



e=0 



We write 

(per m X)(r) m ) =X m (fj m ) = (X^ T (rj m ),X^ T (rj m ),---,X^ T (r) m )f (m > 0) (2.4) 



wc sec 



and call per m X = X m the perturbation tensor field of order m. From Lemma 2.1 
that 

X m (fim) = (X^ T (f IO ),X^ T (f) 1 ),...,X^ T (fj m )f 

and hence the perturbation tensor field per m X = X m has a characteristic: the z-th com- 
ponent depends only on rjo, r/i, • • • , rji, not on any rjj, j > i. 



Lemma 2.2 Let X G TJ(M) and S G T(M), we /iaue 

(X'(u)[S(u)])«(r),) = ((X(«))»)'(^&], i > 0, 
where Si = (S^ T ,S^ T ,---,S^ T ) T . 



(2.5) 



Proof: We first have 



(x^ymsi] = - 



8=0 



X(ui + sj2e k SW) 



k=0 



d_ 

86 



6=0 



X'{ui) [5(^)1+0^). 
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We apply the above equality to the following Taylor series 



d k X( 



U; 



k=0 



! de k 



6 = 



and then get the required equality fl2,5|). jf 

Evidently, (|2j) implies that for X G TJ(M), 5 G T(M) we have 

(per m (^'(«)[5'(«)]))('7m) = (per m X(7? m ))'(^ m )[(per m S)(7? m )]. 



Lemma 2.3 Le£ X G TJ(M). T/ie following equalities hold for any vector field Sm 
(Sq , Sf, • • • , -Sjy) T G T{Mn), where Si, < i < N , are of the same dimension, 



d l X(u N ) 



e=0 



Proof: First from Taylor series 



N 



N 



Qi 



N 



X'(u N )[J2e j Sj], 0<i<N. 

e=0 j=Q 



e i d' l X(u N ) 



we directly obtain 



i=0 



N 



! de l 



+ o{e N ), 



e=0 



(x(«))'(w)[s„] = Ef(^^ 



i=0 



(vn)[S N }+o(e n ) 



e=0 



On the other hand, we have 



(X(u N ))'(fj N )[S N ] 



d_ 

85 



5=0 



N N 

X(u N + 5 £ JSj) = X'(u N ) [E e^i 



i=o 



i=o 



(2.6) 



These two equalities give (|2.6|) again according to Taylor series. The proof is finished. # 



Theorem 2.1 // the operator $ : T(M) — > T(M) is hereditary, then the operator $jv : 
T(M7v) r(Mjv) defined by 



(per N $)(fj N ) = ^n(vn) 
(ijv(w))i 



i,j=0,l,-,JV 



1 d l -^(u N ) 



(i-j)l de^i 



£=0. 



(JV+l)x(AT+l) 



$(770) 



1! 8s 



e=0 



J_ d N $(u N ) 
W. de N 



e=0 



l_ d$>{u N ) 
1! 9e 



£=0 



(2.7) 



is hereditary, too. 
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Proof: Let K N = (K? , K? , • • • , K?) T , S N = (S? , &[,■•-, S?) T G T(M N ), where K u S u < 
i < N, are of the same dimension. It suffices to prove that 



& n {vn)[$nK n ]S n - $ N & N (r)N)[K N ]S N 
-& N (f) N )[® N S N ]K N + $n$'n(w)[Sn]K n = 0. 



(2.8) 



First we easily get the i-th element of the vector field &nKn and the element in the (i, j) 
position of the matrix ^' n {t)n)[Kn\- 



Kj, 0<i<N, 



(& N (fj N )[K N \) 



e=0 



(w)[Kn] 



£ = 0/ 



1 d { - j 



£ = 



N 

(<f>'(u N )£££ k K k }), 0<i,j<N. 



k=0 



Here we use Lemma 2.3 for the calculation of the second equality. Now we can compute 
the i-th element of ib' N (f)N)[$ n Kn]Sn : 



(& N (f) N )[$ N K N ]S N )i 



= E 

3=0 

i 

= E 

3=0 

i 

= E 

3=0 

i 

= E 

3=0 







ai- 


-3 


(i 


-J)! 


de* 


-3 




1 


Qi- 


-3 


(i 


-J)! 


de i 


-3 




1 


Qi- 


-3 


(i 


-J)! 


de i 


-3 




1 


Qi- 


-3 


(i 


-J)! 


de i 


-J 


— ^ 




1 





N k 



£ = 



*'Me^e 



1 d k - l ^{u N ) 



£ = 



K, 



£ = 



1=0 
N 



K, 



£=0 



£ = 



1=0 

N 



£=0 /=0 



Si 



(i-j-iy. de^i- 1 



$'(ujv) ®(un)Ki Sj, 0<i<N, 



£ = 



0<j+l<i 

and the i-th element of &n&n{vn)[Kn]Sn'- 

(*A(w)[#Jv]Siv)i 



= E 

3=0 

i 

= E 



1 &-i$(u N ) 



1 d^^juN) 

fr< (i ^]y. de>-i 



1 

,ok^y-^ 

3 3~k n 

EE 



TV 



£=0 Z=0 
Qj-k-l 



S k 



ototo^- k ~ i y- d ^ k ~ l 



Q'iu^msk 



e=0 
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i i j—k 

EEE 

k=Oj=k 1=0 
i i—k i 

EE E 



i 



-k-l 



(i-j)Kj-k-l)\ de*-j 

1 &-1&(u N ) 



£=0 



dei~ k ~ l 



<f>'(u N )[Ki]S k 



e=0 



(» - JV-U - k - 1)1 de*-3 

I Qi—k—l 



k=0 l=0j=k+l 
i i—k 



e=0 



-fc-Z 



^(^[i^ 



e=0 



k=0 1=0 

E 

0<fe+Z<! 



e=0 



1 



d 



i—k— I 



(i-k-iy. de^- 1 



e=0 



Sk ) 

(^(u N )^(u N )[Ki]S k ), 0<i<N. 



Therefore by the hereditary property of we find that each element in the left side of 
fl2.8|) is zero, which shows (|2.8| ) holds. The proof is completed. # 

Theorem 2.2 If the operator J : T*(M) — » T(M) is Hamiltonian, then the operator 
J N : T*(M N ) -» T{M N ) defined by 



(per N J){fj N ) = Jn(vn) 



i,j=0,l,-,N 



d l+ i~ N J{u N ) 







(i + j-N)\ de i+ i- N 



■%>) 



e=0J (AT+l)x(AT+l) 



J_ dJ(u N ) 
1! 9e 



e=0 



Am) 

is Hamiltonian, too. 



e=0 



J_ d N J(u N ) 



£=0 



(2.9) 



Proof: Let a N = (a%, af, • • • , a^) T , fi N = ($f , 0[,-", Pn) T \ 7iv = (lo > lT>'", 1n) T G 
T*(Mjv), where 7$ < i < iV, are of the same dimension. We only need to prove 

that 

< Sat, J^(i?Ar)[JAr/V]7Ar > +cycle(aAr, /%, 7at) = 0, (2.10) 

because there doesn't exist any problem on linearity and the skew-symmetric property. 
First noting Lemma 2.3, we can compute the element in the (i,j) position of the matrix 
J'n(vn)[JnPn]- 



(J'Avn^JnPn})^ 

l Qi+j-N 



(i+j-N)\ de i+ 'J- N 



e=0 



N 

J'(u n )[£^(JnPn)i 



(i + j-N)\ ds i+ i- N 



1=0 

N 



e=0 



N N 
J'(u N )[^ £ 



d k+l - N j{u N ) 



1=0 k=N-l 



(k + l-N)\ de k + l ~ N 



ftk 



e=0 
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Qi+j 



-N 



{i + j- 


AT)! 




-N 


1 




Qi+j- 


N 


{i + 3- 


AT)! 




-N 


1 




Qi+j- 


N 



N 

e=0 k=0 
N 



N ^k+l-N ak+l-N jf 



UN, 



e=0 



, s __ k (k + l-N)\ de^- 
j'(u N )[J2e N - k (J(uN) +o(e k ))P, 



N 



e=0 



(i + j-N)\ de i+ i~ N 



k=0 
N 



1 



e=0 

Qi+j-N 



j'(u N )[£e N - k J(u N )(3 k 



N 

N -, 



-N 



k=0 
JV-fc ji 



(e N - k j'(u N )[J(u N )p k }) 



e=0 



E 



Qi+j+k-2N 



k=2N-(i+j) 
Therefore we have 



[i + j + k — 2N) ! de i +i+ k - 2N 



(j'(u N )[J(u N )p k }), 0<i,j<N. 



e=0 



< a.N,J' N {riN)[JNl3N]lN > +cycle(a7v,/?Ar,77v) 

I Qi+j+k-2N 



+cyc\e(ai,/3 k ,jj)) = 0. 



(< Oi, J'(u N )[J{uN)Pkhj > 

e=0 



In the last step, we have utilized the Hamiltonian property of J(u). So the required result 
is proved. # 

Completely similar to the above two theorems, we can show the following result. 

Theorem 2.3 If the operator : T(M) — ► T*(M) is symplectic, then the operator @n : 
T(M N ) -► T*(M N ) defined by 



(per N @)(fj N ) = @n(vn) 

(®N(fjN))i 



i,j=0,l,-,N 



1 d N ~ i ~ j e( 



UN 



j_ d N e(u N ) 
ni ai^ 



j_ d N e(u N ) 

1! dF> 



e=0 



£ = 



{N-i-j)\ de N - i -i 

e( % ) 



e=0. 



(AT+l)x(AT+l) 



J_ d@{u N ) 
1! de 



e=0 



(2.11) 



is symplectic, too. 



We mention that when J and G are invertible and J = © 1 , we have 

(per n J)(w) = ((per^e)^)) -1 , 
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which shows that the inverse of the perturbation symplectic operator (pei N @)(r]j^) is 
Hamiltonian and vice versa. Note that the above new operators have a vector t)n of de- 
pendent variables and hence involve q(N+l) dependent variables when u is a g-dimensional 
vector. The above three theorems also provide us a method to generate new Hereditary, 
Hamiltonian or symplectic operators from a known one. This is interesting in the soli- 
ton theory. In particular, we can put forward the following operators by the first order 
perturbation: 



o 

J(r? ) 
J(Vo) J'(rjo)[vi} 



6(770) 



e(%) 




Taking the first order perturbation once more, we can obtain a little more complicated 
operators: 

heSrft) 



(per 1 per 1 $)(r} 3 



dedS 



(per 1 per 1 J)(?73 







J(Vo) 






J(Vo) 

J'(m)[m] 



e=S=0 


■%>) 


J\m)[m 









(per 1 per 1 G)(??3 



d 2 <~>(ri +ern+5ri2 +e5'q 3 ) 



--8=0 



dedS 

&(Vo)[V2] 

Q'(vo)[m] 
0(770) 

Here we have changed two dependent variables while making the second perturbation. Of 
course this kind of perturbation may be done without any stop at finite steps and hence 
the resulting operators are full of various algebraic structures. 






*'(»&)) fol] 

Avo) 
J'(vo)[m] 
J'(vo)[m] 

d 2 J (rio+ern+5ri2+e8ri3) 
dedS 

6(770) 











6'(77o)[7? 2 



6(770) 





e=5=0 

6(770: 







3 Integrable properties of the perturbation equations 



Let us recall the perturbation series (1.12) 



N 



un = Y1 £ l r]i, fj N = (rjo, 7?f . • • • , 77^) T , N > 1. 



i=0 
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For the evolution equation ( |1.11| ) 

u t = K(u), K G T(M), 

we consider the following iV-th order perturbation equation 

UNt = K{un) + o(e ) or u^t = K(un) (mod e ) 

which leads to the following equivalent equation 

1 dK(u N ) 



f) Nt = K N (f) N ), namely r\ it 



, < i < N. 



£=0 



(3.1) 



(3.2) 



In this section, we would like to discuss integrable properties of the perturbation equation 
( |3.2p , which include recursion operators, iT-symmetries (i.e. time independent symme- 
tries), master symmetries, linear representations (Lax representation and zero curvature 
representation) and Hamiltonian structures etc. 

Theorem 3.1 Let K £ T(M). The operator ® N : T(M N ) T(M N ) determined by 
is a recursion operator of the perturbation equation rjNt = Kn(vn) defined by $3.?\ ) when 
<5 : T(M) — » T(M) is a recursion operator of ut = K(u). Therefore the perturbation 
equation r)N t = KN{flN) has a hereditary recursion operator <&n(vn) once u t = K(u) has 
a hereditary recursion operator $(u). 

Proof: Let S N = (S$, Sf , • • • , Sjf) T 6 T(M N ), where S h < i < N, are of the same 
dimension. By Lemma 2.3, we have 



d k $(u N ) 



de k 



(VN)[K 



3 Tfij) 



Qk 



£=0/ 

&(u N )[K(u N ) + o(e")] 



N 

e=0 3=0 



N .. d k &(u N )[K(u N )} 



Qk 

de k 



£=0 



de k 



0<k<N, 



£=0 



and 



1 d i 



N 



E 



i\ de 



K'(u N )[£e k S k 

k=0 

, 0<i<N. 



£ = 



Therefore from the above first equality, we obtain the i-th element of & N (rj]\i)[KN]S]y: 

1 &S&(fi N )[K(u N )]S. 



(& N (r)N)lKN}SN)i = Yl 



3=0 



, < i < N. 



(3.3) 



£=0 
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From the above second equality, we can compute the following two terms: 



{K' N {r) N )[$ N S N ]) 
i ! Qi-k 

tTL(i- k)l de^ 

k=0 x ' 



i—k 



A 1 & 

£*L (i-k)\ dei-k 

k=0 v ' 
i i -. 

EE 

j=0 k=j 
i i i ai-j 

EE 

j=0 k=j 



(K'(u N )[(<S> N S N ) k \) 

£=0 

k 

k'(u N ){y: 

£=0 j=0 

Qi—k 



1 d k ~ j <f>(u N ) 



{i-k)\(k-j)\ 



[k-j)\ de k ~i 



s 4 



£ = 



K'(u N ) 



£ 



£=0 

k-j 



£ = 



(k-j) 



■K'(u N ) 



de k -i 
d k ~^(u N ) 



S, 



£ = 



de k ~i 



Sj 



£ = 



E 

3=0 

i 

E 

3=0 
i 

E 

3=0 



1 &~3 



(i-j)l de l -i 
1 eh* 



K 



1, £ k-3 d k -i${u N ) 



e=0 k=j K J ' 



S; 



e=0 



(i-j)l de { -3 



£=0 



1 d^jK'iuM^ju^Sj}) 



, < i < N; 



£=0 



(3.4) 



* i d i - k ^>(u N ) 



EE 

i=o fc=j 



E 



1 8*-Jff'(tW)[Si] 



£=0 J = 

d l - k ${u N ) 



^{k-j)\ de k ~3 



£ = 



(i-k)\(k-j)\ cV" fc 



d k ^K'{u N )[Sj] 



e=0 



E 

i=o 



1 &-j($(u N )K'(u N )[S j \) 



de k ~i 



< i < N. 



£ = 



£ = 



(3.5) 



Now in virtue of the above three equalities (|3.3j ), fl3.4j ) and fl3.5|) , we easily see that 



at 



-S N + $ , JV (7? iV )[^iv]^ - K' N (r) N )[$ N S N ] + *jvl^(w)[&r] = 0, 



which implies that $jv(^iv) is a recursion operator of ?77vt = Kn- A combination with 
Theorem ^T] gives the proof of the second conclusion. # 

Theorem 3.2 Let K,S £ T(M). There exists a relation between the perturbation vector 
fields Kn and £V 



[K N (fj N ),S N (fj N )] = (K N )'(fj N )[S N (fjN)] ~ (S N )'(fj N )[K N (fj N )] = f N (f) N ), (3.6) 
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where Tjv(?7jv) *s ^ e perturbation vector field of the vector field T = [K,S] E T(M). 
Therefore we have 

(1) if a E T{M) is an n-th order master- symmetry of the equation ut = K{u), then &n = 

is an n-th order master- symmetry of the perturbation 

equation rjm = Kn{t)n); 

(2) the perturbation equation rjm = Kn(J]n) possesses the same symmetry algebra struc- 
ture as the original equation Ut = K(u). 

Proof: By Lemma 2.2, we see for the i-th element that 

(T(u))W = (Jf(u)[5(u)])W(w) - (S'(«)[ir(u)])(%) 

= (K' N (f ]N )[S N }) l - (S'^fjM^KN}),, 0<i<N, 

which shows (|3.6| ) holds. The rest of the proof is obvious. The proof is finished. # 

The relation ( |3.6| ) implies that the perturbation series ( l-12j ) keeps Lie product of 
vector fields invariant. 



Theorem 3.3 Let K E T(M). When the evolution equation ut = K(u) has a Lax rep- 
resentation (L(u))t = [A(u),L(u)] where L, A are two matrix differential operators, the 
N-th order perturbation equation f]Nt = Kn(vn) has the following Lax representation 

(L N (f) N )) t = [A N (fj N ),L N (fj N )], (3.7) 

where the spectral operator Ln and the Lax operator An read as 

(per N L)(r] N ) = L N (fj N ) 

1 d^iLiuN) 



i,j=0,l,-,N 



(i-j)l det-i 



£=0- 



(jV+l)x(JV+l) 



1 dL(ii N ) 
1! de 







£=0 



1 d N L(u N ) 
N\ de" £=0 

(per N A)(f) N ) = A N (fj N ) 
(A N (f] N )) tJ 



1 dL(u N ) 
1! de 



£=0 



(3.8) 



i,j=0,l,—,N 



d :i -iA{u N ) 



(i-j)l de l ~3 



£=0J (JV+l)x(AT+l) 



A(vo) 

1 dA(u N ) 
1! de 



e=0 



A(vo) 



1 d N A(u N ) 
N\ de 1 * 



£ = 



1 dA(ii N ) 
1! de 



£=0 



Aim) 



(3.9) 
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Proof: We first note that 



u Nt = K{u N ) + o(e N ), 



and thus we have 



Qk 



e=0 



((L(u N )) t - [A(u N ),L(u N )}) = 0, < k < N. 



(3.10) 



Let now us compute the elements of the differential operator matrix [An,Ln]- Evidently, 
we know that [An,Ln\ is lower-triangular, that is, 

{[A N ,L N ])ij = Q, 0<i<j<N. 

On the other hand, when < j < i < N, we can compute 



E 

k=j 



1 d i - k A(u N 



(i - k)\ de l - k 



1 



d k -iL[ 



UN 



e=0 



(k-j)l de k ~i 



1 



E 



i - j\ d l k A(uN) 
i-kj 



de 

1 A{u N )L(u N ) 



i—k 



d k -iL{ii N ) 



e=0 



de k ~i 



e=0 



e=0 



e=0 



In the same way, we obtain 

Therefore we have 

([An, L N ])ij 



O^Liu^AiuN) 



e=0 



1 &-i[A(u N ),L(u N )] 



< i,j < N. 



e=0 



(i-JV- 

Now we easily find that ( |3.7| ) is true due to ( 3.1C| ). The proof is finished. # 

Note that the spectral operator Ln and the hereditary recursion operator <&n have 
the same form of matrix. In fact, we can take the hereditary recursion operators as the 
spectral ones. More precisely, Ut = K(u) has a Lax representationcS Q t = [$, K'\ where 
and K' are a recursion operator and the Gateaux derivative operator of K, respectively. 
The following result for the case of zero curvature representation may also be shown. Its 
proof is omitted due to the completely similar deduction. 



Theorem 3.4 Let K G T(M). When the evolution equation ut = K(u) has a zero cur- 
vature representation (U(u))t — (V(u)) x + [U(u),V(u)] = where U and V are two ma- 
trix differential (sometimes multiplication) operators, the N-th order perturbation equation 
fjNt = KN{fjN) has the following zero curvature representation 

(UN(fj N ))t ~ (U N (fj N )) x + [Un(tin), Vn(tin)] = 0, (3.11) 
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where the two matrix differential operators Un and Vjy are of the form 



(per N U){fj N ) = U N (fj N ) 
(UN(w))i 



i,j=0,l,-,N 



1 #-3U(u N ) 



e=OJ (7V+l)x(7V+l) 



U{vo) 



1 dU(ii N ) 
1! de 



e=0 



U(vo) 



1 d N U(u N ) 
M de N 



e=0 



i a[/(ttjy) 

IT de 



e=0 



(3.12) 



(per N V)(fj N ) = V N (fj N ) 

{Vif(m))ij 



i,j=a,i,-,N 



i d^vim) 



(i-j)l de % -i 



e=0. 



(AT+l)x(JV+l) 



1! de 



£=0 



1 a^Vfav) 



£=0 



1 gV^jv) 

1! de 



e=0 



V(vo) 



(3.13) 



Theorem 3.5 Let K G T{M). If the equation ut = K(u) possesses a Hamiltonian struc- 
ture 

ut = K(u) = J(u) / , 
ou 

where J : T*(M) — ► T(M) is a Hamiltonian operator and H E C°°(M) is a Hamiltonian 
function, then the perturbation equation fjm = Kn{vn) also possesses a Hamiltonian 
structure 

SH N (fj N ) 



fjm = K N {fj N ) = Jn(vn)- 



Sf)N 



(3.14) 



where the Hamiltonian operator Jn^n) *s given by $.9Q and the Hamiltonian function 
Hn 6 C°°(Mn) is determined by 



(per N H)(rj N ) = H N (r) N ) 



1 d N H{u N ) 



N\ de N 

The corresponding Poisson bracket has the property 

i d N 



(3.15) 



£ = 



{H m ,H m } j N (fj 



IN j 



N\ de N 



{H l: H 2 }j{u N ), iy 1 , J ff 2 €C°°(M). 



(3.16) 



£=0 
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Furthermore the perturbation equation fj^t = Kn(vn) possesses a multi-Hamiltonian 
structure 

n Nt = K N (n N ) = Jin{vn) = • • • = JmNiVN)- 



5f) N " 5 • 5fj N 

once Ut = K(u) possesses an analogous multi-Hamiltonian structure 

u t = K(u) = Ji(n) — — — = • • • = J m {u)- 



Su 



5u 



Proof: Let 7 = $M- £ T*(M). Then we have 



1 d\J(u N )j(u N ) 



i\ de l 



e=0 



E 



d j j(u N ) 



e=0 



dei 



, < i < N. 



e=0 



Thus we get 



fjNt = K N {f) N ) = Jn(vn)in(vn), 
where the cotangent vector field 77V G T*(Mn) reads as 



(3.17) 



IN 



(vn) = ( w 



1 c^ 7 T ( 



UN 



Nl de N 



_ '(7V-1)! de N ~ l 



1 dj T (u N ) 



e=0 



1! de 



e=0 



,7 T (%)) 
(3.18) 

We hope that the cotangent vector field 7^ is a gradient field. If so, the potential function 
should be the following 



H N (fj 



N) 



< 7iv(A?)jvMjV > dX 



1 " 1 ^ d^(\u N ) 



1 d N 



Nl de N 



< 



de 1 



, n N -i > dX 



e=0 



f 1 

< ~{(Xu N ),u N > dX = — 



1 8 N H( 



UN 



de N 



e=0 



Actually, the cotangent vector field is a gradient field. We can show that 

8H N {f] N ) 



7n(vn) 



Sfj N 



(3.19) 



According to the definition of the variational derivative, we have for any Si G T(M(rji)) 



5 ( 1 d N H(u N 



< 



Srn \N\ de N 

1 d N 



Nl de N 



e=0 

H'iurfpSiirH)] 



,Si{rh) >= ( 



1 d N H(u N ) 



Nl de N 

1 d N ~ i 



(vi)[Si(m)} 



e=0 



(N — i)l de N ~ 



£=0 

H'(u N )[Si(TH)] 



e=0 
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1 d 



N-i 



(N-i)\ de N ~ 



1 d 



N-i 



£=0 



5un 



(N-i)\ de N ~ 



< 7(ujv),Si(?7i) > 



£ = 



< 



1 d N - i j(u N ) 



(N-i)l de N ~ 



,Si(m) >, o < i < n. 



£=0 



This shows that Q3,19| ) holds, indeed. Therefore fl3.17| ) is a Hamiltonian equation. 

Let us now prove the property (|3l6|) . Let = ^,7 = G T*(M). In virtue of 
( |3.19|) , we can compute that 

{H 1N ,H 2 N\j -Avn) =< ,Jn{vn) — > 



Sfj. 



N 



Sfj 



■N 



N 



1 d N - i P{u N ) 



^ < (N-i)\ de N ~ i 

d i+ j- N j(u N ) 



N 

E 



i 



Ai + j-N)\ 



Q e i+j-N 



e =0 j=N-i 

1 d N ~ij(u N ) 



£ = 



(N-j)l de N ~i 



> 



N 



i d N - i p(u N ) 



1 d N 



£=0 

i d l (j(u N ) 7 (u N )) 



N\ de N 



e=0 

< f3(u N ),J(u N )^(u N ) > 



1 d N 



£=0 



N\ de N 



> 

£=0 

{H x ,H 2 }j(u N ). 



£=0 



It follows that the equality ( 3.16| ) is true. 

Further noting the concrete form of new Hamiltonian operators, a multi-Hamiltonian 
structure may readily be established for the perturbation equation. The proof is completed. 

# 

We should note that two important formulas: ( 3.15 ) and ( |3.18| ). ( |3.15| ) provides a 
explicit formula for computing constants of motion of the N-th perturbation equations 



and ( 3.18| ) gives rise to an expression of perturbation cotangent vector fields. 



4 Applications to integrable equations 



4.1 MKdV hierarchy 



We first consider the following MKdV hierarchy^ 



5H 

u tn = K n = a nx = J^ n a Q = J-^ 1 , n>0, 

ou 



with 



J = d, *(«) = --d' z +ud~ 1 ud, H n (u) 



bn+l — Cra+l 

2(2n + 1) 



, n > 0, 



(4.1) 



18 



where aj, 6j q, i > 0, are recursively defined by 

a = n, 6 = 1, c = -1, 
a i+ i = ^aj, j > 0, 

= \a ix + d~ 1 uda i , i > 0, 
c i+ i = ^a^ - d^udai, i > 0, 

and dd^ 1 = d~ 1 d = 1. The first equation is exactly MKdV equation 

1 3 2 



(4.2) 



Its inverse scattering transform was first studied by WadatiE§ . The MKdV hierarchy ( pLl] ) 
possesses zero curvature representations HI' 

U tn -V^ + [uM n) }=0, n>0, 

with 



u = 


u —A 






-1 -it 





n > 0. 



(aA»)+ (6A") + A 
(cA")+ -(aA»)+ 

Here the plus symbol + stands for the choice of non-negative power of A and 

oo oo oo 



(4.3) 
(4.4) 



i=0 



i=0 



i=0 



Applying the integrable theory in Sec. ||, we can obtain infinitely many new heredi- 
tary operators per^per^ • • -per^^, N\, N2 • • • , N m > 1, starting from the hereditary 
operator of MKdV hierarchy 

4 

We easily find that some explicit expressions: 



(pei N $)(fj N ) = diag( ~d 2 , ■ ■ ■ , ~d 2 ) + 



N+l 



fc=0 



(7V+l)x(iV+l) 





drjod 1 r]i + drjid 1 r/o "' 



dr] d 1 7] 1 +drjid 1 7] -jd 2 + dr] d lr q 



(per 1 per 1 <l>)(?7) = 

" -l^ + dmd ~i m 
dr]Qd~ 1 r]i + drjid^rjQ -\d 2 + dr) d~ li i]o 








\d 2 + dr] d 1 T] 









Efc=o%;<9 i r/ 3 _ fc 0ry o «9 L rj 2 + dmd L Vo drj d % + <%d L f] -\d 2 + drj d L r] 
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where r] = (rjo, rjx, ry 2 , r/ 3 ) T . It needs a large amount of calculation if we directly prove the 
hereditariness of the above two operators. 

The 7V-th order perturbation equation of ut n = K n reads as 

flNt n = (pei N K n )(rj N ) = (per Ar $)(77 A r)(per Ar if n _i)(57 A r) = ((peT N $){f) N )) n f) Nx , (4.5) 

among which is the N-th. order perturbation equation of MKdV equation (f4.2|) 



Viti 



i i 



1=0 fc=0 



1 



-rVixxx + ^ X VjxVkm, 0<i<N. 



j + k + l = i 
0<j,k,l<i 



(4.6) 



The resulting local evolution equations defined by (f4.5|) are all integrable soliton ones for 
any n,N > 1. They all have zero curvature representations 

((per N U)(fj N )) tn - ((per N V^)(fj N )) x + [ (pex N U)(f) N ), (per N V^)(fj N ) ] = 0, 

and bi-Hamiltonian formulations 

VNt n = {per N J)(r] N ) = {pei N M)(r] N ) — 



SfjN " ' dVN 

where the Hamiltonian operator M = JL = — jd 3 + dud~ 1 ud. Therefore they possess 
infinitely many common symmetries per^i^m, m > 0, and infinitely many common con- 
stants of motion peT^H m , m > 0. In particular, for the iV-th perturbation equation 
( f4.6[ ) of MKdV equation ( |4.2j ) we can get the following explicit results: two Hamiltonian 
functions 

1 N 



(per N H )(f) N ) = »^2vjVN-j, 



N 



(per iV ii"i)(?7 A r) 



1 



N 



12 



J2 'U'lX .i..r.r I ^yX! VjxVN-j,x 



3=0 



+ - 



E 



3=0 
r lji r lj2 r ]j3 r ]j4.'l 



ji + 32 + 33 +ji = N 
< 31,32,33,34 < N 



a Hamiltonian pair 



(perjvMXf/jv) 



(per N J)(rj N ) = diag( d, d, • • • ,d ] 

N+l 





-\d 3 + drjod 1 rjod 
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and a pair of Lax operators 



(peT N U)(fj N ) 



A 



(per N V<V)(rj N ) 



A N ■■■ Ax Aq 
where the operators Ai, Bi, < i < N, are given by 



Bo 

B\ B 

Bn • • • B\ Bq 



A 



, Bi 



-i -m 

Vi^ ~ jVixx + 2 E ji 4- j 2 +j 3 = i Vjl^Vja 



-8 i0 x + \{r) ix - Y.)=o mn 



<$ioA 2 + 2 (Vix + Ej=0 VjVi-j)^ 
< ji,h,33 < i 



In addition, we can generate a r-symmetry algebraic of the perturbation equation (|4[ 
by a perturbation of r-symmetries of ut n = K n . Moreover we may also consider the 
nonlinear problem of the Lax systems 



4>Nx = {pei N U)(fj N )4> N , 4> Ntn = (per N V^)(f) N )cf> N , 



similar to Ref. 



4.2 Kadomtsev-Petviashvili equation 

Let us now consider the Kadomtsev-Petviashvili equation 

u t = K(u) = K 3 (u) = d~ 1 u yy - u xxx - 6uu x . 

It has Lax pair 

L(u) = —j=d y + d x + u, A(u) = 3id x l Uy — 3u x — 6ud x — Ad x , i = y/^-L 
v3 

and time independent symmetries 



(4.7) 



(4.8) 



K\ = ~u x , K 2 = u y , K n = [K n -i,T], n > 2, 

where r = yK% + \xu v + ^d~ l u y and further the k-th order master symmetry (see for 
example Ref. pF 



T (fc) = yk-i^ r W ^ = ad^ad^ • • • ad Ki .y k+j -\ J, k > 1. 

Therefore KP equation ( |4.7| ) possesses the following time polynomial dependent symme- 
tries 



1=0 
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which contain the symmetries with the forms J^o - ^ and J^lf^...^-- 

The corresponding perturbation equation fjm = becomes 

i 

T) it = d~ l r) iyy - rj ixxx VjVi-j,x, 0<i<N. 

j=0 



(4.9) 



Following the general theory in Sec. y, ( [4.9| ) is also an integrable equation in 2 + 1 dimen- 
sions. For example, it has Lax pair (Z n (r/jv))t = [An{vn), Ln(vn)], where the spectral 
operator Ljy and the Lax operator A^ read as 

m -^d y + <^ + % 



L 



A' 



VN 



m ^9 y + 52 + 7/0 



A 



N 



A(vo) 
Md x 1 r] ly - 3rji x 



3id x 1 7] Ny - 3r] Nx - 6r] N d x ■ ■ ■ 3id x 1 T]\ y - 3r] lx - Qi]id x A(j] ) 
Moreover (^1]) has the k-th. order master symmetry 

per 7V T( fe ) = (/-\0 1 _^) T , fc>l, 

N 

and thus possesses time polynomial dependent symmetries 

k t i 

pei N a {k) = J2 -n(adper iV x) Z per w r (fc) , k > 1, 
1=0 L 

and 

k j.1 

^ 7T (ad p er Jv ^) Z per iV 4 / ; ) 2 ... i j,k> 1. 



(k) 



1=0 l - 



Although KP equation has no regular recursion operatorfU, we may also construct a bi- 
Hamiltonian formulation of ( |4.9| ) through the perturbation of the bi-Hamiltonian onel^3 . 
But here we omitted the discussion because of the complicated notation. 

Let us now take another form of the Kadomtsev-Petviashvili equation 

Uxt = d x K = u yy - (u xxx + Quu x ) x . (4.10) 
Evidently its solutions include all solutions of the Kadomtsev-Petviashvili equation (ET 



and hence it is a little more general than (|4.7|) . We are about to see that there exists a 
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different symmetry property between ( [4.7] ) and Q4.10 ). Now the corresponding linearized 
equation to ( [4,10 ) is as follows 



&xt = d x (K'[a]) = dyyO - d^a - 6dl(ua). 



(4.11) 



One may directly prove that (4.10) possesses the following three time dependent symme- 
tries 

1 



^ (2) (/) 
- (3) (/) 



6 



ft + fu x , 



Tzfuy + ^ftyu x + fu y , 



l 
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= --^fttty 2 + ftt{\y 2 u x 

60 O 

2 2 1 

+ ft(^ u + ^y u v + % xu x) + Mi 



(4.12) 
(4.13) 

(4.14) 



with an arbitrary function / of t, each of which is not any symmetry of fl4.7[ ) while ftt ^ 
0, fttt 7^ or futt 7^ respectively. That kind of symmetries is first introduced in Ref. 
pi]]. In fact, there is a rule to generate these symmetries. For example, 



^ (2) (/) 



^ftty + ft[K, ~y] + f[K, [K, -ly]]. 



We can also construct new vector fields 

n jn-i 



i=0 



dt r 



n > 4. 



(4.15) 



But in these vector fields, there exist nonlinear terms involving u y or u etc., which can't 
be balanced in the linearized equation ( 4.1 1] ) . Therefore ( [4.15| ) are not any symmetries of 
Q4.10| ). The symmetries determined by ( 4.12| ), ( |4.13| ) and ( [4. 14|) constitute a Lie algebra 
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(4.16) 



where /, g are two functions of t. Note that the second commutator relation is somewhat 
different from one given in Ref. [32]. But if we choose a {2 \f) = VZcr^Kf), then two alge- 
bras are the same. The above symmetry algebra contains the following three subalgebras 
{<jW(/), aM(f)}, {a^if)} and{ffW(/), a {2) (f), <r (3) (1) = u t }. The last subalgebra has 
the simple commutator relations 



k,aW(/)] = a«(/ t ), [ut ,a^(f)]=a^(f t ). 
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It is easy to see that the iV-th order perturbation equation of ( f4.10| ) 

i 

Vixt = Viyy ~ {Vixxx + 6 ^ VjVi-j,x)x, <l< N (4.17) 

3=0 

has also the similar symmetry property. More precisely, it possesses the following time 
dependent symmetries with an arbitrary function of t 

perjvcrM (/) = (-lf t + frj 0x Jr] 1 x,---,fVNx) T , 

per iV cr( 2 )(/) = (--^ftty + ^ftyVOx + fvo y , \ftyr}\ x + fr]iy, ■", \ftyr}Nx + fVN y ) T , 
per N aW (/) = (-^fttty 2 + ftt(^y 2 V0x - jgx) + /t(|r/ + jyvoy + \x-q Qx ) + fr] 0t , 
ftt{\y 2 riix - jgx) + ft(lm + Ivviy + l x Vlx) + fvu, 
■ ■ ■ > ftt(ly 2 VNx - jgx) + /t(f »7iV + fywy + I^Wx) + fvm) T , 

which constitute the same infinite dimensional symmetry algebra as ( [1.16 ). In fact, the 
equation ( [4.17| ) has the same integrable property as ( |4.10| ), It is an interesting problem 
how to construct integrable equations which possess a kind of symmetries involving an 
arbitrary function of time variable. 

It is well known that the evolution equation has not a similar property, i.e. it doesn't 
possessf^II the following symmetries involving an arbitrary function / of the time variable 

n Jn—i f 

f-i dt n ~* y h 

i=0 

where Si, < i < n, don't depend explicitly on the time variable. In general, it possesses 
time polynomial dependent symmetries generated by its master symmetries'". Therefore 
Q4.7D and ( 4.10Q have different symmetry algebras. 



5 Concluding remarks 



We may also make another perturbation series 

N e i 



(5.1) 



i=0 



similar to the perturbation series ( |1.12| ). This moment, the corresponding perturbation 
equation reads as 

(5.2) 



where the perturbation vector field is of the form 



,< if _ dK T (u N ) 



Oe 



d N K T [ 
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de N 



e=0 
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In an analogous way, we can generate another new hereditary operator 



$iv(6v) 



dehi 



(5.3) 



e=0J i,j=0,l,---,N 



from a known hereditary operator $(n), and &n(£n) is a recursion operator of ( |5.2| ) pro- 
vided that <&(u) is a recursion operator of the original equation U( = K(u). From the above 
expression, we see that the formation of the perturbation equation under the perturba- 
tion series ( |5.l[) is simpler but the corresponding recursion operator is more complicated 
than ones under (1.12). However there exists an intimate relation between two kinds of 
perturbations because we have a Miura transformation 



m = 7?;(£o,£ir--,6v) = < i < N. 

il 



(5.4) 



For example, the operator &n(£,n) may be generated by a transformation 



1! 



2! 







1! 



2! 
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-i 



(5.5) 



In fact, we can similarly obtain any new tensors by means of the Miura transformation 
(|5.4D - For instance, a new Hamiltonian operator may be engendered by 



Jn(€n) 



dfj N 



Jn(v 



N) 



djN 

dfjN. 



d l+ i- N J(u N ) 



(i + j-N)\ de i+ i~ N 



(5.6) 



OJ i,j=0,l,-,N 



We remark that by the resulting perturbation equations in Sec. ||, we can generate 
approximate solutions of the original equations to a precision o(e N ). This is different from 
the construction of the r functions in bilinear formation, where the expansion series holds 
exactly for any order precision. It is also of interest to note that the perturbation equations 
are all integrable coupling with the original equations and the original ones always appear 
in the first position. Therefore our integrable theory provides an approach for constructing 
integrable coupling of soliton equations and enriches the intention of perturbation bundle 
established in Ref. fjq| . However it is still a problem deserving of investigation how to 
construct more general integrable coupling by perturbation. 
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